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Abstract 



Quadratic theory of gravity is a complicated constraint system. We investigate some 
consequences of treating quadratic terms perturbatively (higher derivative version of backre- 
action effects), which is consistent with the way existence of quadratic terms was originally 
established (radiative loop effects and rcnormalization procedures which induced quadratic 
terms). We show that this approach overcomes some well known problems associated with 
higher derivative theories, i.e., the physical gravitational degree of freedom remains un- 
changed from those of Einstein gravity. 

Using such an approach, we first study the classical cosmology of R+ (5B? theory coupled 
to matter with a characteristic p oc a(t)~" dependence on the scale factor. We show that 
for n > 4 {i.e.,p > |p) and for a particular sign of /?, corresponding to non-tachyon case, 
there is no big bang in the traditional sense. And therefore, a contracting FRW universe 
{k > 0,k — 0,k < 0) will rebounce to an expansion phase without a total gravitational 
collapse. 

We then quantize the corresponding mini-superspace model that resulted from treating 
the PR^ as a perturbation. We conclude that the potential W{a), in the Wheeler DeWitt 
equation — ^ + 2W{a) ip{a) = 0, develops a repulsive barrier near o ~ again for n > 
4 {i.e.,p > ^p) and for the sign of (3 that corresponds to non-tachyon case. 

sign of /?. 

Since a ~ is a classically forbidden region, the probability of finding a universe with 
singularity (a = 0) is exponentially suppressed. Unlike quantum cosmology of Einstein's 
gravity, the formalism has dictated an appropriate boundary (initial) condition. 



Classical and quantum analysis demonstrate that a minimum radius of collapse increases 
for a larger value of \f3\. 

It is also shown that, to first order in (3, I3B? term has no effect during the radiation 
(p = |p) and inflationary {p — —p) era. Therefore, a deSitter phase can be readily generated 
by incorporating a scalar fleld. 
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I. INTRODUCTION 

Since the discovery of the singularity theorem of Hawking and Penrose [1] , the speculation 
of creating non singular theory by incorporating quantum property of gravity and/or using 
non- Einstein gravity has attracted some interest. Because of advancements in quantum field 
theory, the two avenues of speculations seem to be mathematically related. That is, even 
if one starts with Einstein's gravity, renormalization consideration dictates that action for 
gravity must have terms that are quadratic in Ricci tensor [2]. 

In quantum cosmology, canonical quantization is the preferred formalism. This is because, 
for a universe that is homogeneous and isotropic in large scales, there is no asymptotically 
"in" and "out" fields, which is necessary in order to implement covariant quantization for- 
malism. The task of identifying dynamical degrees of freedom for quadratic gravity has been 
reduced to solving constraint system by Boulware [3] . Because of the technical difficulties of 
solving such constraints, quantum cosmology for quadratic gravity has been solved for only 
simple systems like vacuum [4,5]. 

those 

In our paper, we explore consequences of viewing R + gravity as R+ perturbation 
(higher derivative version of backreaction) [6]. In essence, in this approach, the physical 
gravitational degree of freedom is not changed from that of Einstein gravity. It is shown 
that this view overcomes the bulk of technical difficulties with the higher derivative content 
of quadratic gravity. In section (II) , we argue that this view is also in agreement with the way 
the existence of quadratic terms were originally established (via renormalization procedures 
that treats quadratic terms perturbatively) . 

Using such an approach (section III), we study the classical cosmology of (5R? theory 

coupled to matter that has a characteristic p oc a(t)~" dependence on the scale factor. We 
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show that for n > 4 {i.e.,p > |p) and for a particular sign of P, corresponding to non- 
tachyon case, there is no big bang in the traditional sense. And therefore, even for a close 
FRW metric, the universe will rebounce without a complete collapse. 

In section (IV), we quantize the corresponding mini-superspace model, which resulted 
from treating the as a backreaction. We conclude that the potential W{a), in the 

ip{a) — 0, develops a repulsive barrier near o ~ 
again for n > 4 {i.e.,p > |p) and for a particular sign of /5, which corresponding to non- 
tachyon case 

Sign conventions used in this paper are as follows, g — (—,+,+,+), R^^, — ^g/^^i^R — 
(+) SttGT^,, (+) R{fx, u) = V^V, - V,V^ - V[^,,]. 

II. PRELIMINARIES 

The most general quadratic action for gravity coupled to matter is 



Wheeler DeWitt equation -|^ + 2VF(a) 



IGttG 



Imatter + surface term 
(2.1) 

We have formally included a surface term to cancel any boundary term that would result 
in applying the variational principle. By dimensional analysis. Pi, (32,^3 arc dimcnsionless. 
We will be interested in applying the formalism to homogeneous and isotropic metric, i.e., 
Wely tensor vanishes Cabcd = 0[2]. By definition of Weyl tensor, CahcdC"'^'^ = RabcdR"'^'^ - 
2RabR''^ + \R^. This gives one relationship among the possible quadratic terms. 

The second relationship is from the 4 dimensional generalization of Gauss Bonnet formula 

[2], 
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E? — ARahR""^ + RabcdR 



exact derivative. 



(2.2) 



The two relationships, combined with the fact that Euler Lagrange equations are un- 
changed by addition of an exact differential, allow any two of /3i,/32,/53 to be set equal to 
zero in the action (2.1). We choose to set [3^ = 132 = 0. 

Upon variation of the metrics, the resulting Euler Lagrange equations are 



which reduce to a familiar forms for = 0. 

For /3 = case, it is well known that both the left and the right side of (2.3) vanish under 
covariant derivative (the right hand side by local conservation of energy momentum tensor 
and the left hand side by Bianchi identity) . It is interesting to note that this is true even for 
(3^0 (the additional term is also covariant constant by the virtue of Bianchi Identity) . This 
consistency is expected since local conservation of energy momentum tensor and Bianchi 
identity are intimated related to reparametrization invariance of the action (2.1). 

In this paper, we shall consider only the simplest metric, that of spatially homogeneous 
and isotropic universe (i.e., FRW metric) 




(2.3) 



The trace of this equation is 



6 ■ IGtiGPR^ + R = SttGT, 



(2.4) 




(2.5) 



with the standard energy momentum tensor 



Tab=P9ab+{p + p)UaU, [/° = 1, W^O,i 



= 1,2,3. 



(2.6) 
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Using, 



for Ricci scalar, one gets for the time-time component of (2.3). 

\aj 3 \2 a a"* y 3 

and for any one of space-space components, 

2-+ - +^-167rG/3 -i?' + 2i? - + 2 - +— -2 , M = -SttG^. 2.9 

As is evident, for /3 7^ 0, (2.8) is transformed from first order to third order, and (2.9) is 
transformed from second order to fourth order differential equation. 

Even on a classical level, there are pathological problems with higher derivative theories 
[6,7]. One of which is the need for additional initial conditions to completely specify a 
system, and whether solutions obtained by solving (2.8 -2.9) for /3 7^ have well behaved 
properperties as /3 — > 0, and existence of run-away solutions. 

A review of motivation for studying quadratic gravity (i.e., renormalization consideration) 
seems to offer an alternate method of handling higher derivative theories of gravity. 

In perturbative covariant quantization, even if one starts with the Einstein's action cou- 
pled to matter as the bare action. 



'6are 



^ I V^R + Imatter (2-10) 



by quantum loop effects from both gravity and matter, the effective action acquires terms 
with higher derivatives 
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^effective — hare + j OiiR^ + a2-Ra6-R"'' + CUsRabcd.R"'^'''^ (2-11) 

with divergent ai,a2,CK3. The precise nature of divergence depends on a choice for a 
matter field. 

The perturbative renormahzation prescription is to add terms to the bare action to 
precisely cancel these infinities, i.e., 

+ + (2.12) 



'bare 



' countei — term 



= J [(A - + (/32 - a2)RabR'' + (/?3 - a^)RabeciR''"^] ■ (2.13) 



This renders the resulting effective action finite, which may be used for semiclassical 
analysis. (For a fuller discussion of unitarity and renormalizability, please see [2,7].) The 
crucial observation is that physical degree of freedom for asymptotic "in" and "out" fields 
were those of Einstein action. Moreover, even when the bare action had higher derivative 
counter-terms (2.12), the higher derivative terms were treated perturbatively and not on an 
equal footing with the Einstein's action. This is in mark contrast with [3]. The details of 
our proposed procedures will be shown in section III. Our proposal for gravity is similar to 
Jaen's et al.[6] procedure for reducing a general Lagrangian system with arbitrary higher 
derivatives into a second-order differential system. 

A further support of this interpretation of theories in which higher derivative were induced 
can be found in classical electrodynamics. For a pedagogical review, please see Barut[8]. 
Consider a non relativistic harmonic oscillator. It obeys the Newton's law. 



mR = Fext = — ma;^R. 
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(2.14) 



If the particle is also charged, than the accelerated particle emits radiation and in turn 
must effect change on the motion of particle. One can take into account of this radiative 
loss of energy by an effective radiative backreaction force 



mR = Fext + Frad (2.15) 

with 

Frad = (2.16) 

The resulting equation of motion is changed from second to third order. Besides the 
necessity of additional initial conditions, (2.15) also has an unphysical run away solution. 
For example, if Fgxt — 0, than there should be no acceleration and hence no radiative loss, 
i.e., R = 0. Yet, it is straightforward to verify that (2.15) admits an unphysical solution 



R = Ce^ with r 



2e^ 



For a weak radiative loss, there are well known approximate methods of handling these 
problems, which result in physically acceptable solutions. For SHO problem, weak radiative 
loss means that the solution should be harmonic to first order, R ^ — w^R. Upon substi- 
tution, the resulting effective equation of motion is returned to the original second order, 
i.e., 

K + co^rR + io^R^O. (2.17) 

For a general Fext, one can still eliminate unphysical solutions by using an integration 
factor and surgically choosing initial condition for R(0).ln either case, the lesson is that the 
induced higher derivative forces were treated as backreactions which did not increase the 
physical degree of freedom. 



The analogy is even stronger for action at a distant treatment of classical electrody- 
namics. In this case, the backreaction force is not deduced by balancing energy but by a 
dynamical process. Here the backreaction force can be split into contributions from near and 
far field produced by the particle. The backreaction force from the far field is relativistic 
generahzation of (2.16), and backreaction force from the near field results in classical mass 
r enormalization . 

The similarity between induced higher derivative (via radiative processes) in classical 
electro dynamic and the present problem with gravity is obvious. Therefore, using these 
as motivations and possibly even as justifications, we shall also treat the higher derivative 
terms perturbatively when implementing the canonical quantization procedure (i.e.. Wheeler 
DeWitt equation). The rest of the paper can be categorized as consequeces of such an 
approach. 

III. CLASSICAL EVOLUTION OF R + R'^ GRAVITY COUPLED 
TO MATTER 

In (2.3-2.4), we are interested in treating contributions from i?^ as perturbation. We will 

use /3, as a dimensionless expansion parameter and study the first order correction to equation 
of motion. The value of /? is of course unknown, but in order to implement perturbation 
method, we will also have to assume that /? is small. Hopefully, results obtained by assuming 
small (3 will only be amplified by a larger (3. We will return to this issue at the end. 

By (2.4), R = SttGT + 0{p) and = -SttG (t,^ - \Tg,,) + 0{(5) 

Therefore, Eq (2.3) is 

^-Rga^ - Rab = GT^b + (^GT'g,, - 2GTT^ - 2r;^g^ + 2T„.,) + 0{(5''). (3.1) 
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We have introduced a notation G = SnG. 

Using (2.5-2.7), to first order in /3, the time-time component is 

0' + ^ = ^Gp - 2G'(3 (1g(3p -p){p + p) + 2%i3p - p)) . (3.2) 
And any one of the space-space components gives the single equation 

- (S) ^ ^ ^'/^ (^(3. - .) (. ^ .) - ^^Mi^^) . ,3.3) 

The matter sector must satisfy local conservation of energy momentum tensor 

^ = -3pa\ (3.4) 
da 

For P — 0,it is well known that solutions obtained by (3.2) and (3.4) automatically satisfy 
(3.3). As pointed out in section (II), this is guaranteed even for 7^ by reparametrization 
invariance of the action (2.1). Therefore, we shall proceed to solving (3.2) with (3.4). 

Prom the form of (3.2), one can immediately extract several conclusions. First, to first 
order in /9, the R"^ has no contribution in the radiation era because of Prad = —^Prad ■ Second, 
if energy momentum tensor ever becomes dominated by an almost constant potential of a 
scalar field, (— ~ p</) ~ V^(0) ~ constant) then we can again conclude that the contribution 
from i?^ vanishes (to first order in (3). Therefore, a possible deSitter phase in R + I3R^ with 
a scalar field should be more or less identical with a deSitter phase in standard Einstein 
gravity. It is of course a separate question whether one can generate an inflationary phase 
in R + R^ gravity without fine-tuning. 

In comparison, Mijic[9] and Page[10] have studied the large j3 range and concluded that, 
even for a vacuum, gravity alone can generate an inflationary phase in R + R^ gravity. 
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Now let us assume that during any epoch in the evolution of universe, the universe is 
dominated by a matter with a characteristic dependence on the scale factor (i.e., p = ^), 
n=3 for matter era, n=4 for radiation era, and n=0 for an inflationary era, etc. From 
conservation of energy and momentum tensor (3.4), we get p — -^fg-- 

Therefore 

(3p - p) (p + p) = ln{n - A)po' (3.5) 

and 

Imp - P) - -An - 4) (^) ' ^ ^ -„(„ - 4) (- * + \g,) ^ + 0(«. (3.6) 
After simplifying, we get 

+ 2U{a) = (3.7) 



2U{a) ^k(l + An{n - 4)^^/3^) - (l + 3n(n - ^)G'P^) + 0(/3^) (3.8) 

which can be interpreted as an equation describing a particle with a unit mass in a 
potential U{a). For /3 = 0, the form of U{a) is shown in Figure 1. As expected, depending 
on the sign of curvature of three space, (A;>0,A;<0,A; = 0) the universe evolves either as a 
bound state, unbounded state, or as a critically opened state, respectively. 

Now, for /? 7^ 0, the contributions from depend crucially on the sign of /?. First, we 
note that /3 < corresponds to non-tachyon case. This can be readily deduced from (2.4). 
As noted by [3-6], (f) = R evolves like a scalar field with — q.^q^q^ - Therefore, P < 
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is needed to eliminate tachyons. Mijec et al.[9], Stelle[ll], Teyssandier and Tourrenc[12], 
Barrow and Ottewill[13], Mazzitelli and Rodrigues[14] have noted that /? < is necessary 
otherwise the Hubble parameter grows without bound. . 

Confining ourselves to /3 < 0, notice that for n > 4:,{i.e.,p > |p), U{a) develops a 
potential barrier near a ^ Q (Figure 2). The interpretation is straightforward. For such 
a case, a contracting FRW universe {k > 0,k = 0,k < 0) will rebounce to an expansion 
phase without a total gravitational collapse. The graph of region p > |p is shown in Figure 
3. a. For comparison, the strong energy condition for isotropic and homogeneous system [15], 
which predicts an existence of singularity, gives p + 3p > and p + p > 0, Figure S.b. As 
is obvious from Figure 3. a + 3.b, most of the parameter space {p,p) which is predicted to 
have singularity in Einstein gravity do not have singularity in i? + jSR"^ gravity. 

Several comments are in order. First, strong energy condition comes from study of 
Raychaudhuri's equation, which describes how a congruence of timelike geodesices deviate 
from one another. Indeed, the appropriate strong energy conditions for R + R"^ is different 
from that of Einstein gravity and can be shown to be in agreement with the present result 
[16]. 

Second, Page[10] and Coule + Madsen[17] have studied a vacuum model and noted a 
bounce solution for only k < 0. The difference with the present work is more than a vacuum 
versus a nonvacuum model. 

As explained in detail in section II, in our approach, we have treat the /3R^ as a backre- 
action on Einstein gravity [6] . In essence, the physical gravitational degree of freedom has 
not been changed from that of Einstein gravity. This approach has the advantage of having 
a smooth limit as /? ^ 0, and avoids the pathological situation with the necessity for "extra" 
initial conditions in higher derivative theories. 

11 



On the other hand, [10,17] method is straightforward yet the field content, gravitational 
degree of freedom, is different from that of Einstein gravity, and the existence of limit as 
/? ^ is questionable. 



The prediction from classical analysis is indeed interesting, but near the Planck time, a 
quantum analysis is needed. There have been some literature on quantizing mini-superspace 
model for R + R^ gravity [4,5]. In the literature, because of technical difficulties of quantizing 
higher derivative theories, only very simple cases (i.e. no coupling to matter) have been 
considered. With our proposal of treating term as a perturbation, we can readily address 
more realistic systems with matter. 

Even though R? contains terms with higher time derivative than in the Einstein's ac- 
tion, we are treating it as a perturbation. Therefore, the physical degrees of freedom (e.g., 
canonical momenta) are determined by Einstein's action alone [6]. 

For FRW metric, the only gravitational degree of freedom is the scale factor a(t). There- 
fore, canonical momentum conjugate to a(t) is [18] 



IV. QUANTIZATION OF MINIS UP ERS PACE MODEL 



5Ig 



3^3 



(4.1) 



47rG'/fc3/2 



aa — 



aa. 



In terms of TTa, the time-time component of (3.2) is 



TTa^ + 2W{a) = 



(4.2) 




(4.3) 
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U{a) was defined in (3.8). 

We will not need the expression for Wheeler- De Witt equation in its most general form. 
To quantize the system, we replace TTa — > —i-§^ in (4.2) to get 

ip{a) = 0. (4.4) 

Because we are interested in semiclassical analysis, we have neglected the factor ordering 
parameter [19-21]. A comment is in order. Rigorously, there should also a "kinetic" term 
^ representing quantum fluctuation of some matter field. Indeed, the analysis gets quite 
involved and it will be addressed in subsequent work. Here, we shall be satisfied with 
particular "semiclassical" analysis in which only the gravitational sector is treated quantum 
mechanically [18, 22-25] 

As before, we will assume that during any epoch in the evolution of universe, the universe 
is dominated by a single matter with p — For example, scalar field conformally coupled 
to gravity would be n=4, a massive quantum field would be n=3, and during a possible 
inflationary era ~ ¥{(()) ~ constant or n=0, etc. 

For /? = and closed {k > 0), the resulting Schrodinger equation resembles quantum 
mechanical description of unit mass in a bound state potential W{a) [18,22-24] (Figure 4). 
The form of W{a) for a universe which has undergone a standard inflationary phase (via a 
scalar field) is shown in Figure 5. Since a = is boundary of a physically allowed region, a 
boundary condition (initial condition) ■0(a = 0) must be specified to completely describe a 
system [26-29]. 

Now, for /? 7^ 0, the situation is similar to what we have discovered by classical analysis. 
Notice that for again (3 < and for n > 4 {i.e.,p > |p) W{a) develops a potential barrier 
near a ^ (Figure 6). 
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region. 

Physical interpretation of '0 is unclear in quantum cosmology, but if one may interpret 
{ip oc Exp{—\ I)) as an indication of small probability, then our analysis indicated that a 
model with p > |p and /3 < will have a very small probability of a big bang or total 
recoUapse. Notice that unlike Einstein's gravity, there is no issue of boundary condition. 
There is a shortcoming. In Figure 6, the dotted region is precisely when the first order 
approximation breaks down. 

VII. CONCLUSION 

In this paper, we have studied R + [3B? by treating W' term as a perturbation. For FRW 
metric with matter, using classical analysis, we have shown that for n > 4 {i.e.,p > |p) and 
for a particular sign of P that corresponds to non-tachyon case, there is no big bang in the 
traditional sense. And therefore, a recoUapsing FRW closed universe will rebounce without 
a complete collapse. 

From quantization of corresponding minisuperspace model, we have shown that the po- 
tential W{a), in the Wheeler DeWitt equation — ^ + 2W{a) ip{a) = 0, develops a repulsive 
barrier near a ~ again for ?i > 4 {i.e.,p > |p) and for the sign of P that corresponds to 
non-tachyon case. Since a ~ is now strictly a classically forbidden region, a probability of 
finding a universe with singularity (a = 0) is exponentially suppressed. 

And in closing, we can address the effects of a larger value of (3. From (3.7-3.8) -|- Figure 
2 and (4.2-4.3) -|- Figure 6, the minimal classical radius of collapse and the size of classically 
forbidden region increase with larger \(3\, respectively. 
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FIGURE CAPTIONS 

1. U{a) for A; > and k < universe, P — 0.. 

2. U{a) n > 4,/3 < 0. 

3. Graph of region for p > |p, Figure 3. a. Graph of regions p > — |p and p > —p, for 
strong energy condition, Figure 3.b. 

4. W{a) for /c > universe, (3 — 0. 

5. W{a) for very early universe during inflationary era, n ~ 0, (3 — 

6. W{a) forn > 4,/3 < 0. 
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